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^ ■ Abstract 

■ In anticipation of a subsequent application to QCD, we consider the case of QED at high tem- 

perature. We introduce a Fradkin representation into the exact, Schwingerian, functional expres- 
sion of a fermion propagator, as well as a new and relevant version of the Bloch-Nordsieck (BN) 
model, which extracts the soft contributions of every perturbative graph, in contradistinction to 
the assumed separation of energy scales of previous semi-perturbative treatments. Our results are 
applicable to the absorption of a fast particle which enters a heat bath, as well as to the propaga- 
tion of a symmetric pulse within the thermal medium due to the appearance of an instantaneous, 
0^ ■ shock-wave-like source acting in the medium. An exponentially-decreasing time dependence of the 

incident particle's initial momentum combines with a stronger decrease in the particle's energy, es- 

' timated by a sum over all Matsubara frequencies, to model an initial "fireball" , which subsequently 

00 ' 

' decays in a Gaussian fashion. When extended to QCD, qualitative applications could be made to 



Oh: 



X 



RHIC scattering, in which a fireball appears, expands and is damped away. 
PACS numbers: ll.lO.Wx, 11.15.Tk, 12.20.Ds 

Keywords: Finite temperature, Bloch-Nordsieck, soft photon, Fermion damping, QED, Matsubara formal- 
ism, Matsubara frequency, functional method, plasma. 
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I. INTRODUCTION 

In a previous article a hot quantum field toy model was used in order to test and 
appreciate the calculational efficiency of the functional methods long developed and used 
by one of us in a large variety of situations (for example, Refs. 21,0, 4|). The result for 
the two-point function came out both non-trivial and remarkably simple, opening on some 
interesting physical interpretations. The scalar model, however, had little to connect it to 
the physical theories of QED and QCD. In the case under consideration of QED, and at 
the same level of approximation, the same two-point function exhibits a far richer structure 
of entwined contributions and associated mechanisms, which are the matter of the present 
article. 

The mechanisms for depletion of a high-energy particle's energy E and momentum p 
when incident upon a medium at equilibrium temperature T, where p ^ T, suggested in 
the following sections are intended to be a small improvement to the seminal work of Wei- 
don B. TaUshiba fl. and B.ai.ot and lancn ^ Q of a decade ago. Ou. .eeKni,ne. 
and points of view are somewhat different from theirs, but it is essential to begin by ac- 
knowledging our debt to these authors, who first introduced and implemented the idea of 
a Bloch-Nordsieck (BN) approximation in order to estimate the Physics of such energy and 
momentum depletion. 

It may be useful to note the interpretations that we bring to this subject, and we here 
enumerate the special aspects of our approach which the interested reader will encounter 
below. 

1. At the very beginning, we separate and discard (the infinities of) those aspects of free- 
particle mass and wave-function re-normalization from the specific effects of the medium on 
the particle. This is simple to perform in a functional approach, but rather complicated in 
the conventional, momentum-space expansion of the proper self-energy part of the inverse 
fermion propagator (of the particle which is entering the medium). Specific effects of the 
effective mass change in this model due to the motion of the particle in the medium may be 
found in the Ph.D. thesis of one of us |10! |. 

2. We introduce a modified BN approximation appropriate to the case when the particle's 
momentum is decreasing as it moves into the medium. Energy losses are calculated by the 
usual, Matsubara replacement of — > cj^, and thermally averaged according to the Martin- 
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Schwinger/Matsubara formalism we use, as described in detail in a previous publication 
But momentum loss is described by a separate, "Doppler" mechanism, which replaces the 
constant BN momentum p by p{t) = p(0) exp (— Ft), where t is the duration of time the 
particle has been in the medium, and F is specified by a simple, semi-classical argument. 
We feel this choice of semi-classical BN momentum is more physical than the conventional 
procedure, appropriate to high-energy scattering, of retaining the constant value which the 
particle has upon entering the medium; but we place no particular emphasis on this Doppler 
mechanism for calculating p{t). Other models may well be better for the description of this 
decreasing BN parameter, but this one is simple, and physically reasonable, and has the 
interesting consequence of modeling the appearance of a "fireball" at the initial stages of 
the particle's thermal history. 

3. We view the medium as an effective mechanism for the loss of a particle's energy and 
momentum, without requiring the particle to remain continuously on its mass shell (this 
is good Quantum Mechanics, because the experiments we are describing do not measure 
this property!). Only after thermalization, when p{t) has decreased to the order of T, and 
its derived exponential decay law is no longer relevant, only when the particle joins its 
many identical twins in the equilibrium distribution at temperature T, can the particle be 
supposed to be on its mass shell. 

4. We rigorously maintain the nature of our BN approximation, with all real or virtual 
kf^ coupled to the incident particle required to satisfy \p\ ^ \k\. As a result, all integrals 
are finite, and easily approximated. For reasons stated in Section IIIII (after Eqs. (13. lip 
and (13.161) ). we do not employ the conventional HTL analysis to describe pair-production 
generated by virtual photons emitted in the medium by the incident particle; rather, we 
estimate such pair-production using a straight-forward functional representation, and find it 
multiplying the ordinary Bremsstrahlung (a contribution to the decay exponent of {p'^Y) 
by a factor of g^\n. (p^/m^). Were the coupling large, rather than that of QED, this term 
could be suppressed by the unitary denominator factor as described after Eq. (I3.16p . 

5. We are able to provide an explicit expression for the time-dependence of the 
thermalization process, as the particle's [thermal average propagator^ initially increases - 
corresponding to the "fireball" -and then decreases rapidly, as given by a specific, Gaussian 
decay. Were we to restrict the final Matsubara sum to = only, that fall-off would be ex- 
ponential; but we are able to sum over all ra, and the result is a stronger, Gaussian approach 



3 



to thermalization. 

6. Our model calculation is able to distinguish longitudinal and transverse components 
of the "fireball". We do not actually compute distributions which resemble a true fireball; 
rather, we use the word to represent a short-lived enhancement of probability as a function of 
time in the medium, corresponding to the incident particle's ability to generate a longitudinal 
burst of secondary particles and photons. By "transverse fireball" is meant a short-lived 
enhancement of probability as a function of the incident particle's time in the medium, 
which can serve to generate a symmetric pulse of secondaries in an arbitrary direction. 

The paper is organized as follows. Section [III describes the essential features of our BN 
derivation for the fermionic two-point function, in quenched approximation. The result turns 
out to be remarkably simple. For the sake of completeness, the theoretical steps which come 
before that treatment are deferred to Appendix |Xl Then, a " Doppler" model for the fast 
particle momentum damping is used to conclude Section [III and the Doppler model itself 
described in Appendix [Bl In Section IIII[ the approximation of quenching is removed so as 
to take fermionic loop leading effects into account. The transverse fireball will be discussed 
in Section [TVl A summary and a discussion of our results are presented in Section [Vl 



II. QUENCHING WITHIN THE BLOCH-NORDSIECK APPROXIMATION 
SCHEME 

The main steps of the approach SbTQ clS follows, as succinctly as possible. 

Inherent to the BN approximation scheme, ordered exponentials which appear in the 



rigorous Fradkin representation of the fermion propagator 



a So ds' a F{y-u{s'))\ ^ ^^2.1) 



are suppressed, because, as in all eikonal/BN models, they generate terms proportional to 
soft photon momenta, which can be neglected in comparison to the particle's momentum. 
Because of our suppression of conventional mass renormalization, and the understood ap- 
pearence of spinorial wave functions on either side of the final propagator, the fermionic 
propagator at zero temperature will have its (m — ^7 ■ p) factor replaced by 2m, and will 
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read 



{xo,p\Sc\yo,y} = (2.2) 

poo P 







where = + m^. With the same approximations, the free-fermion thermal propagator 
is given by 

Si,(^,Zo) = (2m)- 5^ / dse-'<-'--'^e-'--^\ 

(2n + l)7r , , 

Un = , zo = xq- yo, (2.3) 

T 

and, by direct evaluation, is equal to 

Sth{co,zo)= (2.4) 
(2m) — {[1 - n(o;)]e-'^"o - n(cu)e+'""°} , 

for 2;o > 0, where n{uj) is the Fermi-Dirac distribution function For the thermal 

propagator in the presence of a background A^-field, the corresponding BN approximation 
gives 



poo 

{xo,p\Gfh''[A]\yo,y) = - J2 (2.5) 



^ ^-is{uj^ -uj^) ^-iuj„zo ^-ig ds' p-A{y-2.s'p) 

where = (u;„,p). The thermal 2-point fermion function, in quenched BN approximation, 
is then given by 

{p,n\sT\y,yo) = e"^ ^ ■ (p,n|Gr[A]|y,yo) ■ Zo[^r], (2.6) 

where Zo[zt] is the free partition function, and in the real-time imaginary-temperature for- 
malism being used for the linkage operator with Matsubara sum, one has D^^ = D^'^ + ^D^^, 
with Y)'^" the causal free-photon propagator, and 5D^^ the proper thermal part of D^^, 



j r _ ik-(u-v) 

= ^'^^(^'-^o)im^D,. (2.7) 



2k e^^ — 1 J t^" 
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with k = \k\. In the Coulomb gauge to be used, one has Aq = 0,V ■ A = 0, and 

D,. = S,, (^S,, - ^) 5,.. (2.8) 

In this first, quenched approximation, all complications related to conventional, T = 
0, free-particle mass and wave-function re- normalization are removed by suppressing the 
{"t / 5J~^c'^ sJ~}'P^^^ '^^ complete linkage operator appearing in Eq. fl2.6p . and by 
retaining only the thermal part of it, that is. 



, (2.9) 

A=0 



g~f / SA^T'^^m'^'IaI: 1 / ^-is(uj^-uil) ^-iujnzo ^-ig ds' p-A{yo-2s'uin,y-2s'p) 

r ^ Jo 

n=—oo " 

which gives, 

T ^ Jo 

n=—oo " 

In essence, this term's contribution corresponds to the particle's energy loss due to the 
bremsstrahlung produced under the enhancement of the heat bath's photons; for ease of 
presentation, that bremsstrahlung produced by the slowing particle inside the heat bath, 
will be re-considered in Section IIIII with the DJ^'^-portion of the linkage operation. 
The argument of the last exponential factor of Eq. (12.101) can be written as 



2., ds, ds,- j 2cos[2(.,-.,)fc.;„]|^p ' ) ' 

(2.11) 

where a factor exp (— |A;|/|p|) has been inserted as a simple way of enforcing the BN approx- 
imation. Certainly, other such limitations are possible and the details of the calculation 
will be changed somewhat, but the Physics will be essentially the same. The integrals in 
Eq. (12. lip are well defined, but cannot be carried out exactly. A sensible, approximate 
evaluation is proposed in Appendix C, leading to an overall form of 

-eV/(-)(p^)V, f{-)^^^ (2.12) 
p p 1 + 1 /p 

where ^ combines some numerical factors and where the approximation for f{T/p) is valid 
in the regime \p\ ^ T. In Eq. (12.121) . it is worth noticing that power of s^; had we obtained 



6 



a power of s instead, then we would be lead back to the scalar field situation in which the 
whole set of BN approximations reduced to an exponential temporal damping of the original 
free-field result (cf. Eq. (3.35) in Ref. |l|). As might have been expected in the case of 
QED, the s^-law leads to a more involved behavior that we now must evaluate. 
One is then left with the expression 

a' = eg\f?f{-) (2.13) 
P 



The free-field result, Sth of Eqs. (12. 3p and (12. 4p . can still be used to re-write Eq. (12.131) in 
the form 



• -t-"o poo 

'^'(^i^)' (2m) -Y ds e-^^('^'-"") e-'""^° = e+^'^J^^' Sth{uj, z^), (2.14) 



and, by using the representation: 



e ^a^' = —— / doe ^a^\ 



it is possible to express our BN-approximated result in the remarkable and rather simple 

form of a Gaussian averaged, cu-translated free-field propagator, 

1 r+00 ~ / \ 

S'r(cu,zo) = ^ y dhe-'" Sth (Va;2_2a6,^o). (2.15) 

However simple in principle, an exact integration of Eq. (12.151) remains out of reach. Fortu- 
nately, the following considerations concerning orders of magnitude are helpful in order to 
extract the essence of the result. 

At small enough coupling constant, one can expect to have < 1. In Appendix O, we 
found ^ itself is small and less than 1. Also, starting from |p| ^ T, all the way down to ther- 
malization, where |p| becomes on the order of T, one has clearly -^/J ~ (T/p) / a/1 + T/p < 1. 
And finally, the essential part of Eq. (12.151) is given by the range of |6| < 1. It therefore 
appears sensible to replace Eq. (I2.15P by 

J dbe-''StH{u;il + g^^/fb),zo). (2.16) 

Inserting Eq. (12. 4p . (with the first, overall factor oi l/u left constant for simplicity), inte- 
gration on the parameter b can be performed, with a convergent and damped result of 



.2m r 1 



u-2T{ 



(2.17) 



As noted in Item 5 of the Introduction, had we retained only the n = term of Eq. fl2.13p . 
the result would have had the same structure, but with a slower fall-off in zq. Defining the 
quantity 



q = e--/^+W2-'r)' cos 



u-2T{ 



2uT' 



then, the squared modulus of Eq. (12.171) reads 



4m 



+ q — q cos [uozq) 



(2.18) 



(2.19) 



whose leading contribution, in view of Eq. (12.181) . is the first term of Eq. (12.191) . 

the meaning(s) of which will be discussed in detail in Section IIIII In order to represent the 
thermalization of such an incident particle, one may introduce the ratio of Eq. (12.201) . taken 
at a given zo-value, to its initial value at zq = 0. 

Let R{zQT,p{zQ)/T,a) = R{zq) be this ratio. One has 



exp 



(2.21) 



l + T/p{zo)J' 

where we have resorted to the approximation of ~ p, which agrees with the ordering of 
scales: p ^ T ^ m. Restoring all of the conventional units, one can write 

ZoksT czo ksT 



ZqT 



(2.22) 



h Ac mc^ 

where fc^ is the Boltzmann constant and Ac the Compton wavelength of the traveling par- 
ticle. We assume a reasonable, semi-classical model for the decrease of p = p{zo), of form 
p{0) exp {—Tzo) with T = FDoppier = C,ac/\c {kT/mc^Y, where ^ is a numerical constant and 
Ac = h/ mc, as derived in Appendix [Bl Then, 

pizo) = piO) e-^-, r = V , (2.23) 



Ac V 

and the rise and subsequent fall-off of R{zq) can simply be read off from the expression 

1 



R{zo) = exp 



(*«r 


2 

\ CZo 


y mc^ y 


1 Ac 



2^a 



CZo 



Ac 1 + T/p{zo) 



(2.24) 



If one assumes a specific form for the "linear" density of equilibrium heat bath photons in the 
Doppler computation, e.g., p(z/) = N{u) = [exp (hu/kBT) — 1]^^, that is, the conventional 
Planck distribution, then the constant ^ can be evaluated, and a specific value of zq predicted 
for when the exponential factor of Eq. fl2.24p vanishes, and the fireball starts to decrease. 
One finds, for example. 



Zq = 2^a Wo{ -2r^a e v 

c r V c po 



(2.25) 



as the time after which a possible fireba. 
branch of the Lambert- W-function 
radius of 1/e, one has 



1 starts to decrease. In Eq. f l2.25p . Wo is the principal 
12| . Provided its argument lies within the convergence 



2(-) 
Po 



exp 



rac' 



< e 



(2.26) 



and then, Wq can be replaced by its series expansion, and Eq. fl2.25p may be approximated 

as 



Zq ~ 2^a— I 1 + exp 

c 



+ 



(2.27) 



Note that O^^akBT/mc^) < 1 is a necessary condition for Eq. f l2.26p to be satisfied, and the 
series expansion of Wq to be reliable. In a system of natural units such a.s h = c = ks = 1, 
this is equivalent to 0{aT/m) < 1, and this condition somewhat specifies, and restricts, the 
amount by which T is assumed to be much larger than m (remember the assumed ordering 
of p ^ T ^ M). One has then for Azo an estimation on the order of a/m. 

III. FERMIONIC LOOPS 



In order to include pair-production as a mechanism for the loss of the incident particle's 
energy, it is necessary to include at least the simplest closed-lepton-loop, whose absorptive 
part corresponds to the probability of pair production. Let us approximate L[A], defined in 
Appendix lAf in the simplest way, as LL4] = f / / Af^{x)'K^'^{x — y)Ai,{y), where the gauge 
invariant representation of K^^, reads [2] 



(3.1) 



-k\5, 



9 



with the renormahzed part of Hr, 

2a 



IT 



duu(l — u) In ( 1 + u(l — u 



(3.2) 



Here, m is the mass of the looping fermions and renormahzation has been performed so that 
n_R(0) = 0. It is a real- valued function of fc^, as long as < A;^ + 4m^, but if k^ > A;^ + 4m^, 
it develops an imaginary part given by the discontinuity of the logarithm across the cut, of 



value 2mQ{kQ — + 4m^). The real part of the logarithm, for large fc^, is proportional 
to ln{k^ /m'^), plus additive constants. 

We have now, instead of Eq. (12.61) . the expression 



(p,^|Sj/, \y,yo} = e 



2 J SAu th SA 



- {p,n\Gth[A]\y,yo) 



L[A] 



(3.3) 



A=0 



where Z[ir] is the normalization factor corresponding to the interacting partition function. 

Using the BN approximation specified in Eq. (12.51) . as well as the simplest approximation 
to L[y4], the functional differentiation of Eq. (13.31) can be performed exactly with the help 
of the functional identity 



exp 



exp 



"2 
■ exp 



■ exp 





6 6 ' 



A-K-A-i 
1 



(3.4) 




f-A 



A=0 



th 



1-K D 



th 



--Trln(l-K-DiO 



where = gp^ ds' S{x — (y — s'p)) and the Trace-Log determinantal factor of Eq. (13. 4p 
has no relation to the traveling particle and is absorbed into the partition-function relation 



m 



e 2 



|Trln(l-K-Dtfe) 



-i/3 



■Zo[ 



The resulting thermal propagator has the same form as given in Eqs. 
except that the term exp {2ig'^ 1 1 P ' ^^tk ■ p} is now replaced by 



(3.5) 

(EH) and fl2J0D . 



exp < 2ig'^ 



.1 f\ 
dsi I ds2 

Jo 



(27r) 



p- D 



th 



1 



KD 



th 



■P 



-2i(si—S2)(kou)n-k-p) 



. (3.6) 
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Recalling Df^ = Y)^'^ + 5Df^, Eq. (12.71) . and fc^ Dc(fc) = D, the denominator in the large 
parenthesis of Eq. (13. 6p . 1 — K ■ Dj^, reduces to 1 + n(/c^), whereas the numerator, with 
Dth{k), separates into two distinct parts, of which we consider first the contribution coming 
from Dc(A;) 

1 



exp 



exp 



exp 



+ n(A;2) ■ ^ 
1 



(3.7) 



1 



2 / /-De-/ 



■ exp 



1 + n 



R 



where the factor of [1 + n(0)] ^ renormalizes all 5^^- dependence in the last line, since Z-^ ^ = 
l + n(0) and g]i = Z^,g\ 

We first consider the first exponential factor on the right hand side of Eq. (13. 7p : In a 
T = 0, non-BN calculation where > p is permitted, this term generates the UV divergences 
associated with mass and wave-function renormalization, and those terms should properly 
be discarded as in Section [Tll In a T > context, this term describes the damping of the 
particle's energy due to ordinary bremsstrahlung (in contrast, the remaining part related 
to SDth in Eq- (13.61) describes the damping of the particle's energy as "enhanced" by the 
thermal photon heat-bath in which the particle is slowing down). 

The evaluation of this first term begins with 



exp <^ iig piPj / dsi 



Sl 



ds' 



d^k g-2is'(feo'^n-fc-p) 



(27r) 



k"^ — ie 



(5, 



(3.8) 



and the integration over ko can be carried out by contour integration: 



dh 



-2is'kouJn 



+ ■ 



k I 



[ko - {k - ie)] [ko + {k - ie)] 
e(n)e-2^"''="" + 0(-n)e2^'''=^ 



(3.9) 



Relying again on the approximations used in Section [TTl which amount basically to the 
neglect of the oscillating factors of Eqs. (13.81) and (13.91) . Eq. (13.91) then becomes simply the 
quantity in/k. Inserted into Eq. (13.81) . one gets immediately a contribution of 



exp 



-—a{p ) s 
on 



(3.10) 



11 



That is, the first term in the right hand side of Eq. (13. 7p adds the amount ^a(p^)^ to 
the a^-constant which appears in Eq. (12.131) . an additional damping independent of the 
thermahzed heat bath (at this level of approximation, of course). 

The remaining factor of Eq. (13.70 displays a nice example of the basic unitarity of QED: 
even though we have used the lowest (7^-order approximation to L[y4], in conjunction with our 
BN treatment, a very large cannot produce an overly large effect, for automatic damping 
(in the Hartree-Fock sense) is provided by its denominator. At first, we shall assume a 
weak effect, and accordingly replace that denominator by 1; and then, subsequently, the 
modifications will be noted when the complete denominator is used. But before we proceed 
with that very contribution, an interesting point must be made concerning the contributions 
attached to the (5D^^-piece of the full D^^^ propagator appearing in Eq. (13.61) . 

The contribution attached to 5D^^ in Eq. (13. 6p may be written in a way similar to Eq. 

(1331) 

1 



exp 
exp 



■ exp 



f 



i + n(p 



f 



(3.11) 



i + n(o) 



"27 ^' i + n(o)i + n(o) 



/ 



where the first factor, in the right hand side of Eq. (13. lip , is that part already calculated 
in Section [TTl which leads to the result of Eq. (I2.13p . The second term would correspond 
to a pair production mechanism, enhanced by the thermal heat bath photons. However, 
since n/j(A;^) vanishes at /c^ = 0, this term vanishes because 5Df/i(/c) is proportional to (5(/c^) 
(see Eq. (12. 7p ). Over the relevant fc-integration range, the leading thermal contribution to 
n(/c), that is the so-called Hard Thermal Loop (HTL) polarization tensor, YL^^^jkn, k), can 
be as large and even 

this leading thermal piece of Il^^^{ko, k), is also proportional to k"^ and therefore does not 
contribute either, to the second term on the right hand side of Eq. (13. lip . 

Returning to the second right hand side factor of Eq. (13. 7p in its lowest ^''^-order, 
exp { — (i/2) f f ■ DJIr ■ /}, one needs to evaluate 



exp < Aig 



dsi I ds' 
io 



d'^k g-2js'(fcoaJn-fep) 



(2vr)^ 



fc2 — is 



-.2 (P • kf 

fc2 



(3.12) 
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and focus interest on the contribution coming from the imaginary part of n/j(fc^), given by 
— (2za;/3) B(fc^ — 4m^). Note that the renormahzation prescription of n/{(0) = ensures 
that there is no singularity at k"^ = 0. And, therefore, the /co-integral over the real part of 
n/j(A;^) receives no contribution [lo| . 

To evaluate the absorptive part of Eq. (13.121) . consider first 



Q{ko-Vk^ + 4m^) 
dkon ; ; ^ (3.13) 

[ko — k + te\ [ko + k — te\ 



I , 1 1 



. dko 



2^ J ^k'^+infi ko — k ko + k 
where the upper limit of A cannot be chosen larger than the particle's available energy 
(which it gives to the virtual photon, which then produces the pair). And even though one 
cannot make a mass-shell measurement of the particle as it passes through the medium, 
its energy surely cannot be too far from its mass-shell value, which is essentially cp (until 
thermalization occurs, cp > ksT). Taking A on the order of p, the oscillating factors of Eq. 
(I3.12P are again sufficiently small to be neglected, and what remains is the simple integral 

IF , 1 1 



k 



, , dk, {- - --— ) (3.14) 

p + k ^/k"^ + Am? — k 
In + In 



p-k a/A;2 + Am? + k 



With k = xp, this combination can be approximately reduced to 

(3.15) 



1 

k 



-.2 

In (1 - a;') + In r ^ 



'4m^' 

Then, combining all factors and retaining only the most important ln(p^/m^) dependence, 
one finds for the absorptive contribution of Eq. (13.121) . the amount 

-^(f)'Apr-in(EL\ (3,16) 



3n \An J \m^^ 

The remaining denominator factor of Eq. (13.71) can be taken into account by writing 1 + 
n_R(A;^) = l+g'^{u—iv), and identifying the new absorptive part of nj:j/[H-nK] as —ig'^v/[{l + 
g'^uY + g^v"^]. Here, one has u = 0{ln{k^ /m"^)), which, after integration, translates into a 
denominator factor of ln(p^/m^), thereby removing the ln(p^/m^) factor of Eq. (13.161) . and 
effectively substituting a factor of [ln(p^/m^)]^^. However, other, higher-order corrections 
from the photon polarization may change the result. From this simple photon bubble, in our 

13 



BN approximation and for small coupling, there appears little pair-production enhancement 
of the Bremsstrahlung damping of Eq. fl3.10p . 

In our calculation, pair-production is considered as one of processes of energy depletion 
of the incident particle. The fermion-loop pairs are not considered as thermalized at the 
instant of production, and have no knowledge of the medium without subsequent interaction, 
which is irrelevant to the incident particle's energy loss. In contrast, previous calculations 
[7, y, O] have replaced internal photon lines with effective (resummed) photon propagators 
in the HTL approximation, in which the closed-fermion-loop momenta are assumed to be 
larger than those of soft thermal photons in the construction of the photon polarization 

nnnn 

tensor [15|, [161, [171, [181]. The loop fermions we use are not taken as thermalized; rather, we 
employ the conventional, renormalized, photon polarization tensor, and extract its imaginary 
contribution as that piece of the calculation relevant to pair production. 

IV. TRANSVERSE VS. LONGITUDINAL FIREBALLS 

This Section develops the interpretation of the resulting fermionic two-point function 
which has been introduced in Ref. Collecting all three damping factors, with p ^ T, 
the a^-constant in Eq. (12.131) now becomes 



3n 



3tt ' 37r ' V m'^ 
2 



1 + 



27rT 
p 




(4.1) 



where it is encouraging to recognize a term of 1 + (^^)^, peculiar to rigorous one- loop 



perturbative calculations 



13| . and Eq. (I2.17P can be re- written as 



.2m r 1 -2 
~ 1 2 



<; _e— ^"-£^^0 _ e-^-i^^'"'^^ cos 



u - 2T(- 



a 



zo]}. (4.2) 



The second term of Eq. (14. 2 p describes the disturbance inside the heat bath which is isotropic 
in the medium. Were one to calculate the spatial thermal propagator, by a Fourier transform 
over p at any given time (before thermalization) in the medium, the cosine factor of Eq. 
(14. 2 p would correspond to the appearance of a disturbance propagating with exponential 
phase factors of [i{p- z — Qzq)\ and z + Qzq)]^ with Q the square-bracket constant that 
appears in the cosine's argument of Eq. (14. 2p : and if the dummy variable p is changed to —p 
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in the integration over the second exponential factor, the result suggests the propagation 
of symmetric, "transverse" back-to-back pulses in any arbitrary direction. This is true for 
the free propagator and interacting propagator, and is to be expected of a thermal Green's 
function which not only describes the effects of an incident particle entering the medium, 
but also contains a description of any "tsunami-like" disturbance originating in the medium. 
[One may think here of the emission by the incident particle of a virtual photon with high- 
energy and little momentum, which immediately decays into an electron-positron pair, which 
then comprise and sequentially generate the corresponding transverse fireballs.] 

The factor, e'^^"^ , in the second term of Eq. (14. 2 p comes from the Fermi-Dirac distribution 

-l i£. _|_ 5,^ 

and the combined factor, -e t^jji^ determines the initial relative magnitude of this 
second term compared to the first. The phase is also different from that of the incoming 
fermion by a shift of 2T{^j^)'^. The square modulus of this second term governs what we 
have called the "transverse" fireball, as defined in the Item 6 of the Introduction, 

a 



2uj _af_/ 2 1 



e T 2^' 



u - 2T(- 



^0 • (4.3) 



Similar to the analysis in Section [III for the longitudinal counterpart which now reads 



2 

m 

— exp 



One can see in Eq. (14.31) that the magnitude of the transverse disturbance rises as the 
incoming fermion starts to lose momentum/energy as l/cu^. However, the exponent decrease 
from its initial value is Gaussian and much faster than l/cu^ as time goes on. Then the 
fireball shrinks in all directions once > (or ZqH''^ > {kBT)^^). Hence, a very simple 
prediction arises from this BN-approximated QED calculation; which is a relative increase 
of the transverse disturbance as zq increases, over a duration extent of Azq = 1/T. 

As compared to Eqs. (12.261) and (I2.27P for the longitudinal excitation (and neglecting the 
modification due to —>■ cP), one would therefore have, because of m > aT, the inequality 

(T) (f) 

Azq > Azq for the durations after which transverse (T) and longitudinal (L) excitations 
quickly decay. 



V. SUMMARY 



Here, then, are all the damping factors estimated in a strict Bloch-Nordsieck (BN) frame- 
work, where there are no ultraviolet divergences and no infrared divergences. The evaluations 
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involving the neglect of weakly-oscillating integrands are of course approximate, but quite 
reasonable; and the results are physically correct in the sense that three sources of momen- 
tum and energy loss are included. It is worth noting that the results are also elegant in that 
they involve the complete thermal propagator and provide a continuous time-dependence of 
the process, a "fireball" growth, followed by thermal decay which is Gaussian, rather than 
a simple exponential. We emphasize that, using straight-forward functional methods, we 
sum over relevant contributions of the thermalized photons in the process of calculating the 
rapid thermalization of the incident particle. In contrast to the HTL approach, we do not 
consider fermion-loop lines which define pair-production to be thermalized. 

In this article, the case of a rapid massive fermion entering a thermalized QED plasma 
has been considered. Through the evaluation of the thermal fermion propagator, our focus 
has been to investigate depletion mechanisms that, in the ideal case of an infinite plasma, 
bring the incident fermion down to thermal equilibrium. The same formalism generates the 
probability requirement for the appearance of shock waves developing in the thermalized 
medium Ij. 

In order to go beyond the limitations of pure one-loop perturbative calculations, the 
present estimate is carried out in a non- perturbative way with a BN formalism associated 
with a realistically-decreasing incident particle momentum. Not only do these approxima- 
tions open the road to tractable calculations, but hopefully, they should also be physically 
relevant to the processes under consideration. 

Within a convenient real-time/imaginary-temperature formalism, calculations are first 
carried out with the help of the quenching approximation, and we find a simple and elegant 
expression in terms of the free, non-interacting thermal fermion propagator. Then, the 
quenching approximation is relaxed by taking leading effects of fermion loops into account. 
And remarkably enough, up to the redefinition of a key parameter, the simple form of the 
"quenched result" is preserved. 

For the incident particle, energy depletion mechanisms are taken to be induced by 
bremsstrahlung and pair production, whereas an intuitive, semi- classical Doppler model 
is formulated to account for momentum damping. In contrast to previous calculations a 
complete time-dependent description of the physical processes at play is obtained. In par- 
ticular, the possibility of longitudinal and transverse shock waves is seen to develop in the 
plasma with different amplitudes and phases. Over two different scales of time duration. 
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both excitations start to increase, and then quickly decay with a gaussian law. This rela- 
tively simple QED analysis was motivated by the experimental runs at RHIC; and it will be 
interesting to learn, in a future investigation, if these relatively simple results also hold in 
QCD. 

Acknowledgments 

H.M. Fried would like to thank the J. Schwinger Foundation for a travel grant 
JSF07080000 which contributed to the research and writing of this paper. 

APPENDIX A: ON THE FUNCTIONAL FORMALISM 



The fully dressed Thermal Green's function is 



S' 



th 



^th 



W\ 



(Al) 



where Lt/i[y4] = Trln [1 — ig{'j ■ A) Sth] and the 2)^ operator of the linkage in the configura- 
tion space representation is 



th 



dx / dy 



(A2) 



SA^ivY 

The thermal fermion Green's function Gth[^] is taken in the conventional Matsubara for- 
malism but a Matsubara representation is not needed and is not used for the thermal photon 
propagator of the linkage operator. In the momentum representation of our formalism, the 
thermal photon propagator is separated into two parts. 



D 



(A3) 



with a corresponding splitting of the linkage operator Eq. (1A2I) . The linkage operation can 
accordingly separate into two steps, the causal (T = 0) and thermal (T ^ 0) part, and the 
order of functional operation can be exchanged. For example, 

eLth[A]- 



s: 



th 



.3 



Gth[A] 



(A4) 



In addition to the Bremsstrahlung effect at T > 0, the first linkage operation, with 
exp[2)^^], will produce factors of mass and wave-function renormalization, so that Eq. (1A4p 
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may be written approximately as 



(A5) 



where mass and wave function renormalizations that have nothing to do with the medium 
are included in the fermion Green's functional, Gt/i,ij[A], and closed-fermion-loop functional, 
Lth,_R[^]- For notational simplicity, the renormalization symbol, i?, will be dropped in the 
following, and the mixed representation will hold 



{p,n\S'th\y,yo) 



(A6) 



{p,n\Gth[A]\y^yo) 



where the " quenching" approximation is used; that is, the fermion determinant is suppressed, 
and Eq. flA6p is replaced (with the subscript Q for " Quenched" ) by 



(p,ra|S't/,|y,?/o)Q 
1 



(A7) 



(th) 



Zo[2r] 



(p,n|Gt/,[A]|^,l/o) 



where Zofir] is the free partition function whose relation to Z[ir] is the following. 



One finds 



■ Zo[zr] 



(p,n|S'j/,|y,|/o)Q 

/•oo 

Jo 



X / d[w] exp ^ / dsi / ds2w{si) ■ h ^{si, S2) ■ w{s2) 



X < m — 27 ■ 



JO 

d^k 



p + g 



Jo 



ds'-bth{k)-[w'{s')-2p] e 



-ik-[w{s)—w{s')]+2ik-p{s—s') 



X exp <; ^ I I I ^^^^ik-Ms2)-wis,)]-2:k-p(s,-s,) [^/^^^) _ 2p] ■ Bthik) ■ [W'{S2) 



^0 



(A8) 



where 



h{suS2)= [ ds'e{si~s')e{s2~s'), 
Jo 

h~^{si, S2) = t;—t;—S{si - S2) 

OS2 OSi 



(A9) 
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and 



e ' 2 



iTrln(2/i) 



d[w] 



(AlO) 



•exp|- j dsi J ds2w{si)-h ^{si, S2) ■ w{s2) 

The Bloch-Nordseick (BN) set of approximations is completed by the replacements, which 
relfect the neglect of momentum fluctuation of magnitude less than p, [w{s') — 2s'p] — > —2s'p 
and [w'{s') — 2p] —2p, so as to get eventually the expression 



{p,n\St,^ \y,yo)Q 

(Zo[ir])"^ [(27r)V] "^^^ ^-''^p-y-^nyo) {m - ■ p} 

/•OO 

xi / dse-''^"''+P'^ 
Jo 



(All) 



which, in the main text, is at the level of Eqs. (I2.10p - fl2.1ip . Note that in passing from Eq. 
f lA8l) to flAllD . the huge parenthesis of Eq. flMjl . 



m — 27 



p + g 



d^k 



ds'Iith{k)-[w\s')-2p]e 



-ik-\w{s)—w{s')\+2ik-p(s—s') 



(2vr)^ 7o 

has simply been replaced by |m — ^7 -p}- That is, the medium generated mass-shift, a long 



noticed phenomenon 



19 



20 



2l| . has been withdrawn from our computation, because in this 



article, focus is placed on the depletion mechanisms acting on the incident fermion, and the 
subsequently-generated shock waves inside the thermalized medium. 

APPENDIX B: ON THE DOPPLER MODEL 



We here propose an elementary, one-dimensional derivation which avoids the problematic 
Lorentz transformation of the temperature T 221], of the p{zq) fall-off. Let p{v) be the linear 
density of photons per unit length, at thermal equilibrium. Then p{v) 5v is the number of 
such photons, of energy hv, per unit of time, in the frequency interval 5v. The momentum 
carried by any photon in that distribution is hv/c. What the charged traveling particle sees 
is a Doppler shift of frequencies: for the photons "approaching head on", and for 
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the photons "approaching from the rear", z/ ^ z/_, with 



\ c — V \ C + V 

From the elementary diagram of QED, let rja be the absorption probability of a photon 
by the fermionic line, where rj stands for some numerical constant, a for the fine structure 
constant. 

The number of photons absorbed per unit of time in the frequency interval is thus 
ar]p{i') Su, the same in either front and rear directions. This allows the calculation of the 
momentum change induced by the process, assuming that all the other interactions, with 
the heat bath thermalized photons, average out to zero. One gets 



that is, 

dt 



-J^ — —P Tooppler, (B3) 
Doppler= ^—1 Y I ^ 



and this is the F-constant that appears in Section [TTl Eq. (12.231) . 



APPENDIX C: APPROXIMATION SCHEME IN THE s-INTEGRAL 



The integral in the exponential factor of Eq. (12.111) becomes. 



2 rs psi p poo 

/ / dn{i-c^) / dkk 

271" Jo Jo J Jo 



cos{2ku!ns') cos{2kp(s') _fc 
^" Jo Jo J Jo i 

where ( = cos6. First carrying out both the Si- and s'-integrals, it reduces to 

cos(2sQ(-))" 



Stt 



°° k e~^/P 
dk ■ 







l-cos(2.gW)^^ 



(CI) 



(C2) 



where Q*-^-* = k{pC ± cUn)- The integral over is a bit complicated and can not be carried 
out exactly. To continue evaluation, observe the oscillating factor exp [—is{uj'^ — u'^)] in the 



s-integral of Eq. (12.101) : when s > 



2\-l 



the oscillating factor effectively 



removes any contribution. The arguments of the cosine factors are 



\sQ 



(±)| 



< S„,,axk\pC ± -C SmaxP\pC ± | < ^.^^2^ ^ 
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Since p ^ m and \(\ < 1, 

\pT00n\ \lT^\ |l±^(2n + l)7r||' 
or effectively IsQ^^-*! < 1. Thus, the arguments of the cosine functions are small and can be 



approximated as 



l-cos(2sQW) _ 1(2.QW)2 ^^^^ 



(Q(±))2 "2 (Q(±))2 
Set X = k/p and the fc-integral becomes 



which yields 



dk—rr -=p / dx — -—p^^ 7=P /(— ) 



-s'^Aff fdCil-C) Hdx 



= -es'g'fi-Wr, (C4) 

where all numerical factors have been combined into = Similar approximations are 
made to derive Eqs. f l3.10p and (13.161) . 
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